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It has been argued that the spacetime of our universe can be accurately described by a perturbed 
conformal Newtonian metric, and hence even large density inhomogeneities in a dust universe can 
not change the observables predicted by the homogeneous dust model. In this paper we study a 
spherically symmetric dust model and illustrate conditions under which large spatial variations in 
the expansion rate can invalidate the argument. 

PACS numbers: 98.80.-k, 95.36. +x 



I. INTRODUCTION 



In cosmology, one models the evolution and observables associated with an inhomogeneous universe of density p{x) 
and expansion rate H{x) by employing a Friedmann-Lemaitre- Robertson- Walker (FLRW) homogeneous/isotropic 
model of density p = {p{x)), where (• • •) denotes some suitably defined spatial average. One then assumes that the 
expansion rate and cosmological observables are those obtained in the corresponding FLRW model. 
f^ , One of the proposals to explain "dark energy" calls into question this long-standing (86-year) procedure. The idea 

^^ ■ is that the expansion rate and cosmological observables of a suitably inhomogeneous universe containing only dust, 
^^ ] if analyzed within the framework of a homogeneous model, seems to behave as if the stress tensor also contains a 
[ — . fictitious negative-pressure fluid {i.e., dark energy). 

f^ ' Although this proposal is conservative in the sense that it docs not involve a cosmological constant of incredibly 

OO . small magnitude (presumably rationalized on some anthropic basis), a scalar field of unbelievably small mass, or an 
^^ ' entirely unmotivated modification of general relativity, it is rather revolutionary because it implies that there is no 
dark energy and the expansion of the universe does not actually accelerate (at least, not in the usual sense). 

At present, the idea that the backreaction^ of inhomogeneities accounts for the observational evidence usually 
. _ attributed to dark energy is more of a concept than a predictive model. However, it is generally agreed that if the 
jrt ■ proposal is to be relevant, nonlinearities are required. 

" " ' There have been many criticisms of this approach. One of them [l|, Q] is based on the claim that even in the 

presence of highly nonlinear density perturbations {Sp/ p 3> 1) the metric for our universe can everywhere be written 
as a perturbed conformal Newtonian metric of the form^ 

d.s^ = a^ir) [-(1 + 2^j)dT^ + (1 - 2^)-fijdx'dx^] , (1) 

where dr ~ dt/a is conformal time, 7y is a metric of a three-space of constant curvature, and ijj satisfies the Newtonian 
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^ By the backreaction idea we mean both that inhomogeneities influence cosmological observables (weak backreaction) and that inho- 
mogeneities influence the average/background scale factor (strong backreaction). Weak and strong backreactions are, however, likely 
connected if we assume that averages should be performed along the light cone. 

^ We will refer to a metric written in the form of Eq. JTJ, satisfying the stated conditions, as the perturbed conformal Newtonian metric. 



conditions \ip\ < 1, {dip/dtl^ < a-'^D'ipD.ilj, and (D'-ipDjip)^ < (DW^ip) DiDjip. The covariant derivative with 
the metric jij is denoted by Di. The usual statement is that in the dust case one is allowed to use the perturbed 
conformal Newtonian metric either in the linear regime [i.e., perturbations of every quantity being small) or in the 
weak-field (Newtonian) regime.^ 

The claim is that if the metric can be written in the above form and satisfies the above conditions, even in 
the presence of large inhomogeneities, any cosmological observable will be the same as the cosmological observable 
calculated with ip = 0, i.e., in the homogeneous/isotropic model. This has been described as a "no-go" theorem that 
backreactions can not in principle account for the observations. 

While it is impossible to consider the most general inhomogeneous solutions, there are spherically symmetric 
inhomogeneous dust solutions, which are not perturbations of Einstein-de Sitter, that can be constructed to give 
obscrvablcs similar to ACDM models. These models serve as a counterexample to the no-go argument. In this paper 
we will show why these models can not be described in terms of a conformal Newtonian metric perturbed about a 
spatially flat background, and attempt to understand the implications for the backreaction proposal. 

Indeed, while it may turn out that backreactions are not the answer, we argue that assuming the results of mea- 
surements of the luminosity distance as a function of redshift usually interpreted as an accelerated expansion, the 
metric describing our universe can not be written in the form of a perturbed conformal Newtonian metric where 
a{t) is calculated from the homogeneous dust model. In other words, if the expansion history of the universe is 
well described by the ACDM model, then perturbing about an Einstein-dc Sitter model by the perturbed conformal 
Newtonian metric of Eq. ([T]) is inappropriate, because Einstein-de Sitter would be the wrong background. This is 
because of large peculiar velocities with respect to the background Einstein-de Sitter space. So if inhomogeneities 
are responsible for the observables usually attributed to dark energy, the universe can not be obtained by small per- 
turbations of the Einstein-de Sttter model. In other words, the reason we interpret the observations as evidence for 
dark energy and acceleration of the universe is that we are comparing the observables to observables computed in the 
wrong background. 

As we will discuss, the reason is that the proper meaning of "peculiar" is "after subtraction of a background Hubble 
flow term." We will argue that large peculiar velocities must be present if the backreaction program works, and the 
peculiar velocities are not related to "local" departures from the Hubble flow that would show up as large velocity 
dispersions. 

As an explicit example, consider the argument of Ref. [3|. They propose that the backreaction of inhomogeneities 
in a dust universe modifies the evolution of the effective volume-averaged scale factor and results in an evolution of 
the volume expansion that resembles a ACDM model, rather than the unperturbed spatially flat dust model. If one 
would write the metric for such a perturbed universe in terms of a perturbed conformal Newtonian metric, then one 
would have to use a{t) described by a ACDM model, not the a{t) from an unperturbed spatially flat dust model. If 
one would attempt to express the metric in terms of a perturbed metric with a{t) described by a spatially flat dust 
model, then there would be enormous peculiar velocities in the Hubble flow. 

We explore the issue by considering the spacetime of a particular Lemaitre-Tolman-Bondi spherically symmetric 
model. In the model we study, the observer is centered in a large (Gpc) underdense region in an asymptotically 
Einstcin-dc Sitter universe. This model has large, but not highly nonlinear, density and curvature inhomogeneities. 
We show by explicit calculation that the metric can not be expressed as a perturbed conformal Newtonian metric 
during the entire evolution. We also show that the reason for this can be traced to large peculiar velocities that 
develop with respect to the Einstein-de Sitter background. 

In the next Section we describe the relevant features of the Einstein-de Sitter and Lemaitre-Tolman-Bondi models. 
In Section IIIII we develop the formalism to express the original Lemaitre-Tolman-Bondi metric in the form of a 
perturbed conformal Poisson metric. In Section IIVI we discuss the reason for the numerical results of the previous 
Section. Finally, in Section |V] we present our conclusions. 

II. EINSTEIN-DE SITTER (EDS) AND LEMAITRE-TOLMAN-BONDI (LTB) MODELS 

In this Section we will establish the notation for the homogeneous, isotropic, and spatially-flat cosmological model 
(the Einstein-de Sitter model). We will then briefly introduce the formalism of LTB spherically symmetric cosmological 
models in the synchronous gauge, and turn to a particular LTB model with values of cosmological observables similar 
to the ACDM model. 



^ But already at second order in perturbation theory the spatial part of the metric differs from the temporal part of the metric simply 
from terms quadratic in the peculiar velocities (i.e., an effective anisotropic stress term is in the stress-energy tensor). 



A. The EdS model 

The EdS model describes a spatially-flat, matter-dominated universe. In order to connect with the LTB solution 
of the next section, we can express the line element in terms of the scale factor a{t): 

ds^ = -dt^ + a'^{t)dr^ + a^{t)r^ dil^. (2) 

The Friedmann equation and its solutions are (ao is the present scale factor-in general, the subscript will denote 
the present value of a quantity and k = SttG) 

H^it) = {d/af ^Kpit)/3 = H^iao/af 
a{t) = ao{t/tof/^ 
p{t) -= Po{ao/af. (3) 

At any epoch, the age of the universe is i = 2H^^/3. We also note that the present distance to the horizon is 
dn = 2H„\ 

B. The LTB model 

The LTB model [J, |5|, |6| is based on the assumptions that the system is spherically symmetric with purely radial 
motion and the motion is geodesic without shell crossing (otherwise we could not neglect the pressure). The line 
element in the synchronous gauge can be written in the form 

ds^ = -dt^ + ^^^^dr' + R" (r, t) dn' . (4) 

Here, the prime superscript denotes d/dr. We will denote d/dt by an overdot. The function W{r) is an arbitrary 
function of r. We will also find it useful to define equivalent functions j3{r) and fc(r), related to W{r) by 

W^{r) = 1 + I3{r) = 1 - k{r)r^. (5) 

One can recover the usual Friedmann-Robertson- Walker metric by the replacements 

R{r,t) -^ ra{t) 

k{r) ^ ±1, or 0. (6) 

In spherically symmetric models, in general there are two expansion rates: an angular expansion rate. H±^ = 
R{r,t)/R{r,t), and a radial expansion rate, Hr = R'{r,t)/R'{r,t). 

Assuming a dust equation of state, the Einstein equations may be expressed as 



i?2 RR' 



^ ' ^ P2 PR/ ^ 



These are the generalization of the Friedmann equations for a homogeneous/isotropic universe to a spherically sym- 
metric inhomogeneous universe. 

Adding the equations results in an equation that can be easily integrated to yield the dynamical equation for -R(r, t): 



^(r,t)^Jpir)^^y (8) 



where a{r) is a function of r related to the density p{r,t) by 



-^(-'^) = i^^M)^-(M) - ^'^ 



Equation ([8]) can be differentiated to yield the dynamical equation for R'{r,t): 

f3'{r) + a'{r)/R{r, t) - aR'{r, t)/R\r, t) 



R'{r,t) = 



2R{r, t) 



(10) 



Since we will eventually express the LTB solution using a perturbed conformal Poisson metric, it is convenient to 
use instead of time a new independent variable related to the EdS scale factor: x = a/ao. We also introduce new 
dependent variables Y = R/ar and Z = R' /a. Therefore, our variables are 



ao 

Y{x,r) = — 
ar 

Z{x,r) = — . 



(11) 



The solution Y — Z — 1 corresponds to the EdS model. The dynamical equations for Y and Z are 



dY_ 
dx 

dZ 
dx 



Y 

X 

Z 

X 



Bix Ai 
y2 ' y3 



1 XB2+A2/Y - AiZ/Y'^ 
2^ ^Bix + Ai/Y 



(12) 



(13) 



where the functions Ai(r), A2{r), Bi{r), and B2{r) are related to a{r), j3{r) and their derivatives by 

a{T)^HlayA,{r) 
a\r)^HlayA2{r) 
P{r)=HlayB,{r) 
P'{r) = HXrB2{r). 

Here, Hq = Hj_(to,r -^ 00) = Hr{to, r -^ cxd) = 100/ikm s~^Mpc~^ is the present EdS expansion rate. 

Since each shell labeled by r evolves independently, we can choose initial conditions to be specified at different times 
for different shells. However, for simplicity we choose to specify initial conditions at the same time for all shells. We 
will choose the initial time to be recombination for the EdS model, x~^ = 1100. The initial value of R (or Y) can be 
chosen freely with a suitable rescaling of r. We choose Y = Z = 1 at x^^ = 1100 to match to the EdS solution as 
r — > 00. 

A particular LTB model is generated by a choice for a{r) and f3{r). For a discussion of the free functions in LTB 
models, see Ref. m. For other works based on LTB solutions aiming at showing the effect of inhomogencities, see 
0, i, i, M, [O, MM, M, M, M, M, M, [H, H III . Our goal here is to study a LTB model with observables similar to 
the standard ACDM model. We will adopt the particular model of Alnes, Amarzguioui, and Gr0n (AAG) [23. l23l. l24l| . 
but with slightly different parameters. The model is EdS at large r, the observer is in the center of an underdense 
(compared to the EdS density) region, and a there is region of density larger than the EdS value to compensate the 
underdense region. 

In the AAG model a{r) and (3{r) are given by 



a(r) 
/3(r) 



H^or' 



tt2 2 2 

Hgaor 



1- 

Aa 



Aa 
~2~ 



1 — tanh 



2Ar 



1 — tanh 



2Ar 



(14) 



Alnes, Amarzguioui, and Gr0n choose h = 0.51, Aa = 0.9, aprp = 1.35 Gpc, Ar = OArp. With that choice of 
parameters, at the origin at present h± ^ hr = 0.65. They found that this choice of parameters reproduces the 
luminosity distance-redshift diagram expected for a ACDM model. That choice of parameters results in a "weak" 
singularity at the origin [l2|, as can be see by the fact that p'{r = 0) 7^ 0.^ While this is not a serious issue (the 



* What Ref. [12| | actually claims is that if there is no weak singularity at the center, then for an observer at the center, dA{m — M)/dz 
cannot be positive at 2 = 0. However, it is possible to have dA(m — M)/dz > for off-center observers. 
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FIG. 1: Upper figure: The density compared to tlie density of tfie EdS model as a function oi v = Hoaor for various values of 
X = a/ao- Middle figure: The functions k — —I3(v)/v'^ and —j3{v) = 1 — W'^[v) as a function of u = Hoaor. Notice that W'^ is 
always very close to unity (an approximation we will employ later). Bottom figure: The current proper distance as a function 
of the comoving coordinate v = Hoaor for the LTB model and the EdS model {Hq^ — 5878Mpc). 



observer could be slightly displaced from the origin or a slightly different density profile could be used) we will employ 
different parameters than AAG to mitigate the weak singularity. We do this simply by the choice Ar = 0.15ro. This 
results in the present expansion rate at the origin of h± = hj. = 0.67. 

The density compared to the EdS density is shown in Fig. [1] Notice that at the origin p' vanishes. The curvature 
function fc(r) is also shown in Fig. [TJ The fact that as r approaches 0, k is negative and p is much smaller than 
the EdS value, implies that the evolution deep in the void will resemble a negatively-curved, empty universe. For 
sufficiently large r, k approaches and p approaches the EdS value, so far from the void the universe evolves as an 



EdS model universe.^ The AAG model can actually be thought of as a swiss-cheese model of holes of radius v^^s- 
(See Ref. [3 for a discussion of swiss-cheese models.) 

In Fig.[l]we display the functions in term of the comoving coordinate label r multiplied by HoaQ. That combination 
is related to the present comoving proper distance in units of Hq , 

TJorfp (today) = i7o / „., . dn ^ HoGq / — — dn ^ — — — dvi, 15) 

Jo W[ri) Jq W(ri) Jq W{vi) 

where v = Hoaor. The spatial curvature term W{r) is close to unity, and Z is of order unity, so to a reasonable 
approximation iJorfp(today) = v, the EdS result. (The choice h = 0.51 implies Hf^^ = 5878 Mpc.) In Fig. [1] we show 
the present proper distance as a function of w = Hf^a^r. 

The evolution of Y{x) and Z{x) is shown in Fig. [51 The initial conditions specified at cc^^ = 1100 are Y = Z ~ 1, 
with an initial underdensity in the interior. Since the underdense region expands more rapidly than average, a typical 
LTB evolution will lead to shell crossing in a rather small dynamical time. The occurrence of shell crossing means that 
the synchronous coordinate system becomes singular and the evolution in that system can no longer be calculated. 
AGG delay the onset of shell crossing by endowing the shells with a large initial infall velocity (negative dY/dx and 
dZ/dx). Thus, it takes the system a while for the underdense regions with large expansion rate to reverse the infall 
and this essentially delays the eventual shell crossing.^ 

Alnes, Amarazguioui, and Gr0n also calculate the luminosity distance as a function of redshift, assuming that the 
observer is in the center of the underdense region (r = 0). The photon geodesic equation for the position of the photon 
as a function of time, t{r) or x{r), is found from the LTB metric |22l |: 

di ^ R'{r,i) dx ^ ^^/^ Z{x,v) 

dr ^l + /3(r) '^ dv ^ W{v) ' ^ ''> 

where we have used the new comoving radial coordinate v = Hoaor. The redshift of the photon, z{r) or z{v), is 

MM 

dz R'{r,i) dz l + zZ{x,v) / x dZ{x,v)\ 

d^^^ '^ ^' ^1 + f3{r) ^ d^ ^ x^^ W{v) \ '^ Z{x,v) dx )' ^ ^ 

Equations (|16p and (|17p are solved with initial conditions z = v = Oatx~l. The luminosity distance is then simply 

dL{z) = (1 + z)'^R{r, i) ^ HodLiz) - (1 + zfxvY{x, v). (18) 

A graph of the luminosity distance as a function of redshift is shown in Fig. [3l The striking feature in Fig. [3] is 
that the LTB model resembles a dark-energy model. AAG show that their LTB model is a better fit to the "gold" 
data set of Riess et al. [25 than the concordance ACDM model (at the expense of additional parameters) . With our 
different choice for Ar, dL{z) in our model is slightly different than the AGG choice of parameters, but nevertheless 
demonstrates a shape characteristic of ACDM models. Note that as z ^ 0, ^^(z) follows the curve expected for the 
empty model. This is because the center of the void is, in fact, evolving as an empty universe at late times. 

Whether our (or AAG's) LTB model is a better fit than ACDM in light of new data is inconsequential. Our goal is 
not to champion any particular LTB model as a cosmological model that describes our universe, but rather to explore 
the failure of the no-go argument that inhomogeneities can not describe (ii,(z) because the spacetime of our universe 
can be expressed in the perturbed Newtonian form. 

We note that in the LTB model there is no dark energy, and since there is only a pressureless component to the 
stress-tensor, no individual fluid element undergoes local acceleration. 

We emphasize that observations can not directly establish that there is dark energy, or even that the universe 
accelerates in the usual sense. They only tell us that high-redshift supernovae are about one-half magnitude fainter than 
expected in the Einstein-de Sitter model. Stronger statements have additional (always unstated) model assumptions! 

In the next section we will discuss the formulation of the model using a perturbed conformal Poisson metric. But 
before turning to the main issue of this paper, we note that although no fluid element in the LTB model undergoes 
local acceleration, one can easily show that in the backreaction prop osal a suitably volume-averaged definition of the 
scale factor has negative pressure (see Refs. [a. KLSl l26l. WA 28 ._ 29l. |30| for a discussion). This is related to the averaging 



issue in cosmology as emphasized by Buchert and Ellis (sll, [33 | 






^ We will refer to the value of r or u when the solution is very near the EdS solution as rg^S or ''EdS ■ 

^ As far as the purpose of this model is concerned (having observables similar to the ACDM model), the synchronous gauge does not 

have problems, and, in particular, the LTB solution is an exactly solvable model with nonlinear inhomogeneities for an extended time 

evolution. 
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FIG. 2: Upper figures: The functions Y and Z as a function of v = Hoaor for the indicated values of a; = a/ao. At the initial 
time of the evolution, x~^ — 1100 and Y — Z = 1. For the other values of x, Y is always smaller than unity and asymptotically 
approaches Y — 1 (the EdS solution) for large values of v. For small values of v the function Z is smaller than unity, then for 
larger values of v it becomes larger than unity, and eventually Z approaches unity from above. Lower figures: The functions 
dY/dx and dZ/dx as a function of v for the indicated values of x = a/ao. The function dY/dx never changes sign, while the 
function dZ/dx does. Note that for x~^ = 1100, the values of dY/dx and dZ/dx have been multiplied by lO"'^. 



III. LEMAITRE-TOLMAN-BONDI MODELS AS A PERTURBED CONFORMAL POISSON METRIC 

The LTB model of the previous section is an example of a nonlinear cosmological model. It is also a model that 
exhibits an observable that is very different from the unperturbed model (EdS). Our goal now is to investigate where 
and why the "no-go" argument discussed in the Introduction fails. 

To get insights, we will express the LTB metric, Eq. (01), in the form of a perturbed conformal Poisson metric with 
small potentials. We will consider a metric perturbed about the EdS (spatially flat) background. In the metric we 




FIG. 3: The luminosity distance as a function of redshift. The curve labeled LTB is the LTB model of this paper, the model 
labeled ACDM is the standard ACDM model with JIa = 0.7, the model labeled EdS is the EdS model, and the model labeled 
Empty is an empty universe model. Rather than Hodh^z), we show the usual difference in the distance modulus compared to 
the empty model. 



adopt, vector and tensor perturbations are absent: the perturbations are given by the potentials "0 and 0,^ 

ds^ = a^{T) [-(1 + 2i')dT'^ + (1 - 2(j))dr'^ + (1 - iip^dn^] , 



(19) 



where r is confornial time, adr — dt. 

When considering a dust model (or, in general, in the absence of anisotropic stress), one usually takes the same 
scalar perturbation in both the spatial and temporal part of the metric, that is, ^ = (f>. Also, it is usually assumed 
that the potentials are time independent. 

As we will soon show, however, we will have to allow the two potentials to differ and be time dependent: this is 
caused directly by the nonlinear structure of the LTB model. This can be seen mathematically by the fact that in 
the case ip{r) = (/>(r), the system of equations will be overconstrained if we are beyond linear order. The cause of why 
ip ^ (p a-nd they arc time dependent is actually one of the main issues we will try to understand in this paper: what 
we want to show is that it is not possible to describe the AAG model by means of the perturbed conformal Newtonian 
metric of Eq. ^ as far as the linearized theory is concerned. 

Of course, this behavior is due to the fact that the AAG model has an initial density and curvature perturbation 
which is not within the linear regime of perturbation theory. Even though this is the specific cause of the particular 
behavior of the AAG model, we will try to learn the general physical cause of that behavior. The AAG model will be 
just a tool to understand the general lesson about the importance of voids not dynamically restricted which, as we 
will explain, are a key-feature to evade the no-go argument discussed in the Introduction. 



^ The perturbed conformal Poisson metric is closely related to the perturbed conformal Newtonian metric of Eq. |[TJ. The differences are 
that in the perturbed conformal Poisson metric the spatial perturbation (</>) need not equal the temporal perturbation (tp), and the 
perturbed conformal Newtonian metric also has the additional restrictions on the metric coefficients discussed after Eq. JTJ. 



We point out that the calculations we are going to show are consistent with the standard result that at linear order 
the potentials are time-independent and identical: ip{r) = ^(r). It is indeed possible to show (see Ref. [201) ^^sd if 
there is no time dependence, then the equations will not overconstrain the scalar potentials in the case ij) = <j). 

In the next subsection we will set up the gauge transformation needed to calculate the potentials ^ and (p of the 
metric of Eq. p9)) in the perturbed conformal Poisson metric corresponding to the LTB metric of Eq. ([4]) in the 
synchronous gauge. 

A. Gauge transformation 

To be able to compare the metrics of Eqs. ((4]) and ([19]), we will express the latter in the synchronous gauge. We 
will use the gauge transformations of the linearized theory: the consistency of this will be discussed later. 
The generic form of the perturbed EdS model, both in Cartesian and spherical coordinates, is [33] 



ds^ = a^{T){-{l + 2i')dT'^ -2^^uJdx''dT+[{l-2(|))5^J+D^Jx]dx'dx^] 

= a^{T) -(1 + 2Vj)dr2 - 2uj'drdT + (l - 2(/> + ^T j dr^ + (l - 20 - ^T ) r^d^"^ 



(20) 



where Dij = didj — -j^i jV^ and T = x" ^ x'/^- Because of spherical symmetry we are including only longitudinal 
degrees of freedom in the metric perturbations. 

Under a gauge transformation the metric perturbations transform according to 

^ = V^-O.C-— C (21) 

a 

Cj = u; + (^ + drf3 (22) 

^ = <^-V/3+^C (23) 

3 a 

X = X + 2/3, (24) 

where ( and f3 express the gauge freedom. In radial coordinates, W'^/S = j3" + 2/3'/r. This transformation corre- 
sponds to a diffeomorphism which keeps the perturbations in Eq. ([20)) small: this is actually the definition of gauge 
transformations in the linearized theory. The new coordinates are given by 

T = T + C 

x = a?— V/3 

f = r-/3', (25) 

where again only the longitudinal component has been taken into account.^ 

In the following, tp, Cj, (f>, x, f, x, and r will refer to the quantities in the LTB synchronous gauge, while the corre- 
sponding quantities without the tilde will refer to quantities in the perturbed conformal Poisson metric. 

We can now choose the gauge transformations in order to end up with ip = uj = 0, that is, to end up in the 
synchronous gauge. Taking in account that we arc coming from a metric of the form Eq. (|19p . which has w = and 
ip^O,we find from Eqs. ^ and ^ that 

i, - 0^^ = a,C+— C=(aC)- (26) 

a 

Co = ^ C = -drP = -a/3. (27) 

These expressions can be expressed in integral form as 

If', . . a 



a{t) Jt a{t) 

-^y_ V(r,ii)dii-3C(r)-(^l-^j-f/3(r), (29) 



Only 13' will appear in the equations since the displacement is dr = — /?'. 



10 

where the overbar denotes the quantity evaluated at t = t , the initial time of the evolution (t corresponds to x~^ ~ 
1100). 

The time-independent integration functions ^(j') and /3(r) express the residual gauge freedom, which we are going 
to discuss shortly. But before doing that, let us make use of the remaining two gauge transformations, Eqs. (|23p and 

dMl). 

Comparing the LTB metric, Eq. ([4]), with the form of the perturbed metric, Eq. ([20]) . we see 

1-2^+^ = ^ (30) 

l-^^-it^-^. (31) 



Now compare the perturbed conformal Poisson metric, Eq. p9|) . with the form of the perturbed metric, Eq. (|20p . 
with result that T = x" — x' /'>' = 0. The general solution to this equation is x = Cir^ + C2, where Ci and C2 are 
constant in r (but possibly time-dependent). In our case we are describing a model where % — > as r — > cx), i.e., it 
approaches the unperturbed EdS solution at large r, so we must choose Ci = C2 = x = 0. Therefore T = 0. 
Since we know (f> and T, we can write 

l-2, + 2/,"-l£c = ^ (32, 

l-2, + 2^-lf<.4l^ (33, 

Since we know ( and j3 in terms of the potentials and initial data for ( and /3 from Eqs. (pS)) and (P5)) . we finally find 
the dynamical equations for the potentials ijj{r,t) and (j){r^t): 



^5^y_ ^"(^,ti)dii+6C"(r)-(^l-^j+2/3"(r)--y_ ^{rM)dt^ 
i^li^r)^ ^"(-'^) * (34) 



These are the dynamical equations to solve for (f>{r,t) and ip{r,t). 

Before solving, it is illustrative to evaluate Eqs. ((34)) and (|35|) at the initial time t = i to find the integration 
functions ((r) and P'{r). The result is 



2 a 



W'^in) - 1 



dri 



(36) 



The fact that C(co) = was used in Eq. ([36]) : all the gauge quantities have to go to zero as r ^ 00 where the 
spacetime is exact Einstein-de Sitter. 

The residual gauge freedom expressed by C,{r) and /3'(r) is therefore fixed by the value of (j) at initial time, that is, 
by the initial conditions of the hole in the synchronous gauge. 

Note that /?' must vanish at the center of the hole because of spherical symmetry. Only (^ is allowed to be nonzero 
at r =: 0. 
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B. Potentials 



With some manipulation we can express Eqs. (|34p and (|35p in a more manageable form. The difference of Eqs. ([34]) 
and (j35p yields one equation, while combining Eq. (|34p with Eq. (|35p multiplied by r and then differentiated with 
respect to r yields a second equation: 



! "Hh) Jl \ r I ' a\ t'" 



r \ T 



cp'{r,t) + ^^l_'yj'{r,h)dh + ^^C{r) = S{r,t), (39) 



where the functions V(r,t) and S{r,t) are defined as 

J_ f R^{r,t) R'\r,t) \ 

^ ' ' 2r V a" (t)^' a2 {t)W^ (r) / ^ ^ 

o. ,^ _ 1 ( R\r,t) R'\r,t) 2R{r,t)R'{r,t) \ 

^'' 2r\a^{tya^t)W^{r) a^t)r J' ^ ' 

Using the fact that as r — > oo, the potentials and their derivatives must vanish, we can integrate Eqs. (j38p and 
with respect to r to obtain 



/ df /■ ^ f / P^'^\ - - 1"°° r°° 

^y_ V'(r-,t2)di2 + 3-(^l-^jC»-/3M = y dnnj V{r,,t)dr, (42) 

V'(r,ti)dii + --C(r) = - / 5(ri,t)dri. (43) 



2 1 /■* , , . 2a- 



3 t Jt 3 i 

Final manipulations involving multiple differentiations with respect to t and integration by parts yield the results: 

1 P°° 
i;{r,t) = -J {rj~r^)Q{n,t)dn (44) 

-I ■ /"OO /"OO 

Hr,t) ^ --- {rl~r')Q{n,t)dn- S{n,t)dru (45) 

Z d J J. J J. 

where Q(r,t) = a2(i)F(r,t). 

Note that the integration functions C{i^) ^^nd P'{r) have disappeared from the final solution. They are, however, 
determined by the value of (l>{r,t) at initial time as seen from Eqs. (|36p and ([57)1 . For completeness we give the 
expressions for C(r) and (3'{r): 

1 Z""" 
C(r,t) = f-r-— y a2(r2-r2)T/(ri,f)dri (46) 

/OO 
F(ri,t)dri. (47) 

Now we are in a position to express the potentials, the displacements in time. At ~ aC(''i Oj ^'^^ radial coordinate, 
Ar = P'{r,t), in terms of the new radial coordinate v = HoaoT. In terms of the variables Y and Z defined in Sec. 
IIIB| the expressions are 

2^P{v,x) = f {vl-v^)M{x,vi)— (48) 

Jv Vi 



Jv «1 Jv Vi 



,At 3 f°^, 2 2N,r/ ^dVl 



(49) 
2=^ = i / {vi^v')N{x,v,):^ (50) 

2^ = rL(x,.r)^. (51) 
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The functions M{x,v), N{x,v), P{x,v), and L{x,v) are 



dY{x,v)Y X fdZ{x,v)Y Z^{x,v) ( A2{v) Ai{v) 



M{x,v) = xr-^^^\ -__ ^-L^ +_^^^[ ' 2 



dx J W^(y) \ dx J 2xW^(v) \Y'^{x,v)Z{x,v) Y^{x,v) 



2x \y^(x,v) 
'{x,v) Z 
dx W'^{v) dx 



,r/ \ ^rr .dY(x,v) Z(x,v) dZ(x,v) , „, 

N{x,v) = Y{x,v)^^^ '! \: (53) 



P{x, v) = Y\x, v) + ^^^ - 2y(x, v)Z{x, v) (54) 

Lix,v) =y'i-^-)~^^, (55) 

where, as before, x ~ a/ao. The evolution equations for Y and Z are given in Eqs. (J12p . The initial conditions are 
r(u,S) = Z{v,x) =: 1 at i~i = 1100. It is clear that for an EdS solution {W{v) = Y{v,x) = Z{v,x) = 1, Ai{v) = 1, 
and ^2(^0 = 3) the sources M[x,v), N(x,v), P{x,v), and L{x,v) vanish. Since we will never allow shell crossing in 
the evolution, each shell evolves separately. 

In the next sections we will present the results for the potentials and coordinate displacements. 

C. Results 

First, we show in Fig. |4]the values of the potentials ipi^i ''^) ^-nd 4>{x, v) for various values of a; = ao/a. There are 
several striking features. First of all, at the initial time, x^^ = 1100 (as well as for x^^ = 500), (p and ip are greater 
than unity for v < 0.3. Also striking is the fact that near v = 0, (j> j^ i/j- For x^^ < 100, (p and tp are smaller than 
unity, but they have different signs. 

In the next section we will explain the reason for this behavior. Here, we note that when the potentials are 
greater than unity, they are constant within the horizon. If we take the initial epoch to be a;~^ — 1100, then vh 
is found by solving Eq. (fT6|) . The value of vh for the LTB metric is not too much different than the corresponding 
value for the EdS model, vh = 2x^1'^. In the EdS model, v^ = 0.06, 0.09, 0.14, 0.20, 0.6, 1.15, and 2 at x'^ = 
1100, 500, 200, 100, 10, 3, and 1 respectively. In Fig. [4]we indicate values of u^ for x^^ ~ 500 and 200. As seen in 
Fig. 21 when the potentials are larger than unity, within the scale of the horizon, they are very nearly constant. We 
will return to this point in Sect. FVl 

Now we turn to Fig. [5] where we show the coordinate differences between the LTB metric and the perturbed 
conformal Poisson metric. First, let us examine Af. At the initial point in the evolution Af/i is much larger than 
unity and positive. So even at the very beginning of the calculation we can not express the LTB metric in terms 
of a perturbed conformal Newtonian metric! This occurs even though the metric is everywhere regular, and density 
perturbations are not highly nonlinear. 

Eventually, when x drops below about x^^ — 100, At/t becomes smaller than unity and continues to decrease. In 
the next section we will discuss the cause of this behavior. 

Now let us turn to Ar. At the initial point in the evolution, x^^ — 1100, Ar is positive. This means that r > f, or 
the coordinate label in the perturbed conformal Poisson metric is greater than the coordinate label in the synchronous 
gauge. For other values of x, Ar is negative (i.e., r < f, or the coordinate label in the perturbed conformal Poisson 
metric is less than the coordinate label in the synchronous gauge). 

Finally, we turn to the issue of peculiar velocities. We first remark that the dust LTB metric of Eq. ^ is comoving: 
there are no peculiar velocities. Again denoting a quantity in the synchronous frame by a tilde, the four-velocity 
is u^ = (1,0). However in the Poisson frame there are peculiar velocities (see also (SJI)- At linear order, under a 
coordinate transformation .t^ ^ i^ — e^, we have (quantities in the Poisson frame do not carry a tilde): 

ui'ix) = u^ix) ~ d^e^u^ix) + e^d^n^ix) (56) 

Uf^{x) = Uf^{x) +df_,e''u^{x) + e'^daU^{x). (57) 

From this we find (f is a unit vector in the radial direction) 

u^'ix) = {l-aC-dC,l3'f)^{l-iP,l3'f) (58) 

u^{x) = {-l-aC + d(:,a^P'f)^{-l-ij,a^p'r). (59) 
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FIG. 4: The gauge potentials %Ij{x,v) and (j){x,v) as function of v for the indicated values of x. The potential tp{x,v) is 
indicated by solid curves and the potential 4>(x, v) is indicated by dashed curves. In the lower half of the figure the values of x 
are indicated only for (j){x,v); the values of x for ■tp{x,v) are, from top to bottom, x"^ = 100, 10, 3, and 1. In the upper half of 
the figure, the values of v at the horizon are indicated by a vertical dotted line for x^^ = 500 (leftmost line) and x^^ = 200. 



As expected, the change in the four-velocities is of the same order as the gauge transformation. At hnear order 



-1 and 



5/^1/^ 



The value of u^ would be interpreted as a peculiar velocity. Now we express the four- velocity as u^ = (7, 7W*). The 
three- velocity is w = \/hrr «'" r — av'^f at linear order. Using this. 7U = a'yv^f — a/S'f, so 



jv — aP'r — vx 



1/2 



N — f. 

V 



(60) 



The result for ^v as a function of v = H^aar is shown in Fig. [6] for several values of x. The large initial values of the 
velocity are the key to understanding why the perturbative description breaks down. 
We now turn to the discussion. 
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FIG. 5: The time displacement Af/t and coordinate displacement Ar/r as a function of v for the indicated values of x. The 
time displacement At is negative for x~^ = 1100, 500, and 100. It is positive for x~^ = 10, 3, and 1. For x~^ = 1100, Ar is 
positive, for all other values of x, Ar is negative. Clearly, Ar = rAr/r vanishes in the limit r ^ as it must. 



IV. DISCUSSION 



The linear gauge transformation we have performed in the previous Section is a useful tool to understand where 
and why the no-go argument discussed in the Introduction fails. In this Section, first we will introduce the general 
physical cause of the breakdown of the description by means of the perturbed conformal Poisson metric. Then we 
will show how this general feature shows up in the model here examined. Finally we will conclude the Section with 
some comments about other models developed in the literature. 
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x-i = 500 - 



x-i = 1100 - 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 



v=Hoaor 



FIG. 6: The peculiar velocity in the Poisson frame as a function of v = Hoaor for the indicated values of x. A negative 
(positive) value corresponds to a radial velocity towards (away from) the center. 



A. Lack of a global background 



The description by means of the perturbed conformal Poisson metric breaks down because of the large AH in the 
AAG model. As can be seen from Fig. [71 



AH I 



-ffEdS 



HEdS 



(61) 



is of order unity. A large AH is a typical feature of models with inhomogeneities, and from this point of view we 
think the LTB model is a useful tool to study the consequences of the evolution of structures in the real universe. 
In particular, the AAG version allows us to follow the evolution of the inhomogeneities for a rather large time span, 
thanks to the large initial infall velocity of the shells (negative value of AH/H^ds for x~^ = 1100, as can be seen 
from Fig. [T]). 

A different expansion rate means a differeirt FLRW solution, and even if the initial AH is small, the two solutions 
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FIG. 7: The expansion rates as a function of v for the indicated values of x. 



(void and EdS) will eventually follow very different evolutions. What we have in mind is the comparison of a 
background FLRW solution in the underdensity (the void region) and in the FLRW solution at large v. We will 
denote the value of v beyond which the solution is approximately EdS as tiEds- We make the point that these two 
patches have different backgrounds, even though the "average" background is the EdS solution. The AAG model 
can indeed be thought of as an exact Swiss-cheese model matched to the density and the (zero) curvature of the EdS 
solution at a sufficiently large radius {v > WEds)- We can imagine populating the universe with other holes of radius 
V ~ I'EdS- On scales larger than any spherical hole, the AAG model evolves as the EdS model. From this p ersp ective, 
in the framework developed by Buchert [S^, the backreaction of the inhomogeneities is zero. (See Refs. [3,[la,[l3| for 
a discussion.) Here, we find a lack of a global background in spite of the absence of a volume backreaction. 

It can happen that at some specific time the background of the void patch is close to the EdS solution. This will 
be, however, a special moment because of the different evolution of the two patches. It will be more typical to have 
a large departure from the EdS flow, which will be manifest in a large AH or in large "peculiar" velocities. These 
"peculiar" velocities are the result of comparing the cosmic evolution of two different FLRW solutions, or, in other 
words, the results of perturbing about the wrong background. 

Summarizing, an inhomogeneous universe featuring voids cannot be studied by means of a global homogeneous 
background, and this holds even if the inhomogeneities do not affect, as it happens in a swiss-cheese model, the global 
evolution. 

B. Discussion of the results 

What we have just discussed is clearly shown by the results obtained in the previous section. To understand them 
it is useful to focus on two key quantities. 

First, it is possible to show (see the Appendix,) that the value of at the origin has the following time dependence: 



, ,'avit) 



(62) 
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FIG. 8: The functions (V f3) evaluated at u = (dashed curve) and (^o = <?!>(« = 0, x) (sohd curve) as a function of a; = ao/a 



where a^ = r^^ ^lr=o describes the scale factor of the void region. In turn, this implies that 0o oc AH±. This 
expression makes clear that (j)o can be thought as the time integral of AH. 

A similar relation holds for f3' (see the Appendix p): a P' :^ R— ar. This implies that a u^ ^ a f3' ~ RAH±, where 
we note that au^, the velocity of the comoving matter in the Poisson frame shown in Fig. [51 is, as expected, almost 
the same as the "peculiar" velocity of the comoving matter with respect to the EdS background: 



R' 



Vpec = / dri— {Hr - -ffEds) 

Jo ^ 



RAH± =RHi 



EdS 



AHi_ 

HEdS 



R H^dS , 



(63) 



where the last relation holds because, as we have seen, AH/H is of order unity. (We have also used the fact that 
W{r) is very close to unity.) 

Expressing these key-quantities as a function of AH, which is the source of the effects, will help us in understanding 
the results of the previous section. 

First, at early times {x~^ < 500) "peculiar" velocities are larger than c, as can be seen in Fig. [6l The reason can 
be understood from Eq. (p5)) : at early times -ffEds is indeed large. This highlights the usefulness of the AAG model, 
which allows us to encompass both early times with a large expansion rate and the present time. This also illustrates 
the interesting fact that even with smaller AH we could have large peculiar velocities at early times when one might 
expect the system to be accurately described as a perturbation of an EdS universe. These large velocities cause the 
potentials to be greater than unity and the perturbed conformal Poisson metric to be unable to describe the metric 
in the void region. In the linear gauge transformation this shows up as the fact that peculiar velocities are of the 
same order as the gauge transformation. This is a general breakdown which affects the perturbed conformal Poisson 
metric even if we perform a nonlinear gauge transformation [36| . 

At later times, velocities become order 10~^c as can be seen from Fig. [51 this is caused by the slowing down of the 
expansion. However, both /3' = Ar and (f> do not depend on AiJ, but rather on its time integral, or, in other words, 
on the different size of the void region compared to the EdS region, that is, on R — ar. This dependence appears in (j) 
because it is the spatial perturbation that has to connect these two differently expanding spatial regions. The reason 
for the breakdown of the linear transformation is, therefore, the growth of P' = Ar, that is, the departure of the two 
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metrics. This shows up directly in (j). The basic problem is that the linear gauge transformation neglects terms hke 
(V^/3)^ compared to terms Hke 4>: see Eq. ([25)1 . As you can see from Fig. [51 this approximation is no longer valid for 
a;~^ < 400. So even if the velocities are small, the linear transformation breaks down because their time-integrated 
effect is large. As a consequence we find, for example, 4> j^ ip. It is however possible to describe the late evolution 
of the model in the perturbed conformal Newtonian metric by means of a nonlinear gauge transformation whose 
potential follows nonlinear equations of motion [36| . 

C. Comments 

We will comment now on the models of Refs. [S^ and [33] • According to Ref. [33] there is no breakdown of 
the perturbed conformal Newtonian metric during the dynamical evolution of the particular (it is important to not 
generalize) LTB model the authors have considered. The model of Ref. [3^ does not deal with the possibility of using 
the perturbed conformal Newtonian metric to describe the spacetime of the universe, but it will nevertheless help us 
to make our point because it is a model which can be described with the perturbed conformal Newtonian metric if 
we perform a nonlinear gauge transformation. 

We start noting that in both Ref. [33] and Ref. [3a] the validity of the perturbed conformal Newtonian metric is 
assumed as a starting point. These particular models are not allowed to have large peculiar velocities. The point is 
that this is a particular choice, and, in our opinion, not general as demonstrated by the dependence in Eq. (|63p of the 
peculiar velocities on AH. We would like to stress again that these "background" peculiar velocities have nothing to 
do with anything that can be measured as a local effect. Therefore, the models of Refs. [33] and [38| over-restrict the 
dynamics of the inhomogeneities, or, equivalently, demand the existence of a global background. For example, in the 
model of Ref. [33], the initial AH of Eq. ([55)1 is exactly zero, and in the model of Ref. [S^] it is AH/ HecLS < 10"^. 
Therefore, these models serve as an illustration to show the circular reasoning at the basis of the no-go argument: if 
we demand that the spacetime metric of the inhomogeneous universe has a global background without any departure 
due to differently evolving regions, then we will find the validity of the perturbed conformal Newtonian metric. This 
is the main point of this work, and the linear gauge transformation is a useful tool to understand it. 

It is interesting to note that the model of Ref. [3^ can confront the SNe data, but its luminosity distance does 
not show a big departure in shape from the EdS one, while the AAG model features a luminosity distance similar 
to the ACDM model. This is mainly due to the size of the hole, roughly 10 times smaller than the AAG one. We 
point out that if we want to have a ACDM-likc luminosity distance in a setup with only one hole and the observer 
at the center, then we will likely have large background peculiar velocities. This might, however, be a typical, but 
not a general feature. To illustrate this we have calculated the luminosity distance for the model of Ref. [S^, but 
with a hole roughly 10 times larger: we found that dA{m — M)/dz < always. This can be seen, again, as due to a 
restriction on the dynamics: allowing larger peculiar velocities we find voids with more general dynamics. We stress 
again the usefulness of the AAG model in showing these general features. 

V. CONCLUSIONS 

The LTB model analyzed in this paper is a very useful toy model to explore the role of large inhomogeneities in 
the determination of cosmological observables such as the luminosity distance as a function of redshift. We employed 
the toy model to demonstrate where and why the no-go argument discussed in the Introduction fails. To this end 
we have performed a linear gauge transformation in order to understand the problems we have to face if we want to 
express the LTB model in the form of a perturbed conformal Newtonian metric, Eq. ([T]). 

We would like to stress that we regard the AAG model as a tool useful to show a general feature needed to evade 
the no-go argument: we focused on the general physical cause of the nonlinear behavior of the AAG model more than 
on its particular initial conditions. 

We came to the conclusion that the way the no-go argument is proved is by assuming, as a starting point, the 
validity of the argument itself. The way to escape the no-go argument is by allowing large "background peculiar 
velocities," where these velocities have nothing to do with, e.g., velocity dispersions in clusters or anything that can 
be measured as a local effect. If indeed inhomogeneities are responsible for the "apparent" acceleration of the universe, 
then the evolution of the expansion rate must resemble a ACDM model. If that is true, then there will naturally be 
"peculiar" velocities if one describes the evolution in terms of perturbations of an EdS model. One cannot simply 
assume peculiar velocities are small and then use that as a basis to show inhomogeneities can have no effect. In the 
calculation of peculiar velocities it is necessary to specify a background solution about which to calculate velocities. 
The peculiar velocities can be small in the ACDM background but not in the EdS background. The velocities appear 
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as a departure from EdS of the Hubble flow, preventing the metric from being written in the conformal Newtonian 
form. 

In other words, to evade the no-go argument we need to free the inhomogeneous models of dynamical restrictions 
(small background peculiar velocities) which impose the existence of a global background. 

Another way to see how the no-go argument of Eq. ([T]) fails is through the equations of motion: in order to have a 
model, e.g., an LTB model, explain the observed departure of dL{z) from the EdS result requires nonlinear equations 
of motion, and this is because of the large background peculiar velocities which can generally occur. By this, we mean 
that in order to end up with a large spatial variation of the background H requires a nonlinear evolution like the one 
of the LTB model discussed in this paper. And this can occur even if it is possible, as in the AAG model, to describe 
an inhomogeneous universe by means of the perturbed conformal Newtonian metric from (let's say) z = 2 to present 
time. This is shown here by explicit calculation. 

Therefore, one must reconsider very carefully the statement that inhomogeneities can not have a significant effect 
on cosmological observables such as the luminosity distance as a function of redshift if arguments are based on the 
assumption that the metric describing our universe can be written in the perturbed conformal Newtonian form of 
Eq. (P). 

Now we enumerate some of our findings. 

1. In the model we consider, even at late times when the density contrast is linear and the velocities are small, 
one can not express the LTB metric in the perturbed conformal Newtonian form of Eq. ([1]) by means of a linear 
gauge transformation. The breakdown of the linear transformation clearly shows the importance of background 
velocities. 

The linear gauge transformation breaks down mainly because of the large Ar = /3' . This is caused by the 
different expansion rates (Aif ) of the void and the EdS solution, as can be seen in Fig. [3 which shows the 
expansion rates as a function of v for several values of x. The result is that the two metrics do not remain close 
enough for a linear transformation to connect them. As discussed in the previous section, the basic problem is 
that the linear gauge transformation neglects terms like (V^/3)^ compared to terms like (j). Summarizing, the 
cause of the break down of the linear gauge transformation is the large Aif , and the underlying cause of that 
is the nonlinearity of the equations of motion. The failure described above can be stated as the fact that the 
linearized equations of motion work only for a limited range of times. 

2. It is possible to find a nonlinear gauge transformation to overcome this problem with a perturbed conformal 
Newtonian metric description of the LTB metric whose potential follows nonlinear equations [36| . These nonlin- 
ear equations (which are actually the exact Einstein equations) at some time will typically give large potentials 
and make the perturbed conformal Newtonian gauge description inconsistent (large At = aC in the LTB model 
examined) because at some point the peculiar velocities will be larger than c if we do not restrict the dynamics, 
as shown by the AAG model. 

3. Although at the present time, in the model we consider, the magnitude of the potentials are small, they are not 
equal at the level of the linearized theory, viz. ^ ^ (f). This could naively be interpreted as some anisotropic 
stress or something else funny with the stress tensor. This is exactly the point that has been made in, e.g., 
Ref. Q: Analyzing an inhomogeneous model in the framework of a homogeneous model introduces spurious 
components and behavior to the stress tensor {possibly even a negative-pressure fluid!). 

4. If we start with a perturbed EdS model, then the growth of perturbations will lead to the formation of large 
voids, filamentary structure, and other nonlinearities. Perhaps the LTB models are useful toy models to study 
the effect of that phenomenon. The model explored in this paper is a good example because it is capable of 
dealing with nonlinearities within a large timespan. 

5. Finally, another point worth mentioning is the physical nature of super-Hubble modes of (f) and ip. At the initial 
time, within each Hubble radius patch, and ip ^^^ very nearly constant, but possibly larger than one (with a 
value of about 20 near the origin). That does not mean that we can simply subtract different values, because 
we are only allowed to do this once. If we do it at x^^ = 1100, then at x^^ = 1 we would have and ip near 
the origin of —20. 

Summarizing, inhomogeneities, in particular voids not dynamically restricted, undermine the existence of a global 
background. 
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The interesting question for future consideration is if backreactions^ of cosmological perturbations can substantially 
distort derived conclusions about the energy-momentum tensor of the universe without allowing dramatic violations 
of statistical homogeneity. The present work addresses the key issue that backreaction effects cannot be assumed in 
principle to be small on the basis of the no-go argument explained in the Introduction and therefore opens a wide 
range of possibilities to tackle the study of the inhomogeneities in the universe. 
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APPENDIX: DETAILS OF THE APPROXIMATIONS USED IN SECT. [IV] 



In this Appendix we will discuss the approximations used in Sect. IIVI Although we have checked numerically that 
they hold, we will present here an analytic explanation. 



1. Approximations for (j>o 



We first want to motivate that the time dependence of (j> near the origin can be estimated as 0o ^ ^{ay/a — 1)^, 

where Uy = r^^ R\r=o- ^^ ^^^'^ ^^ ^^- P^]) . 

The second term dominates the expression for <j>{r, t) in Eqs. (|l5)) or (P^ . so in the limit that W{r) = 1, (f>o{t) can 
be written as 



Mt) 



1 / R i?'\^ dn 



n 



2 \ari 



R 

ari 



n 



dri, 



(A.l) 



where we defined rEdS as the coordinate radius at which the metric is the EdS one (it will disappear from the final 
expression). To proceed further, define a function c{r,t) as 
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(A.2) 



Now define ^1^=0 — (^vf where a^ gives the scale factor for the FLRW solution at the center. Thanks to Birkhoff's 
theorem, outer shells do not influence the center, which can be indeed thought as a homogeneous patch. Expanding 
R{r, t) in the void region as^" 



R{r,t) 



i^,{t) + ay2{t)r'' + 0{r^), 



(A.3) 



we find 



Mt)--i- 

V a 



(A.4) 



^ We do not really distinguish between weak or strong backreaction because we think that what matters are the observables. Assuming, 
indeed, that an inhomogeneous model gives a ACDM-like luminosity distance, then it is not straightforward to understand which 
category of backreaction (weak or strong) the model belongs to. The strong backreaction indeed is about a change of background, but 
this is also the case for the weak backreaction because, as in the case of this example, it "defines" an effective background, the ACDM 
model. Moreover, the change of background of the strong backreaction should be related to observable data on the light cone as it is 
for the weak backreaction. 
^^ In the expansion a„i = 0, thanks to the fact that p'(0, t) = as discussed in Ref. |12|| . 
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where, after integrating by parts, we have made the approximation that the void region occupies almost all of the 
volume of the hole of radius rg^s ■ 

Now we wish to show that c{r,t) depends only weakly on time. To this end, we expand {R/ar) and \{Rlar) ] in 
r around rEds to third order, that is up to the fourth derivative of R with respect to r, and then we plug the resulting 
series in the integrals. The result is 

cr,t)~- + - + — -r^]{r-r^as)- (A.5) 

The first remark is that the term proportional to R"" (rEds^t) cancels. The dependence on time is therefore in the 
R"'{rEdS,t)/R"{rEdSit) term. Following a procedure similar to this one used to obtain Eq. (3.5) of Ref. [11], but 
keeping both density and curvature inhomogeneities, we find: 

R"'{rEdS,t) ^J_^ d"'{rEds) ,^ g. 

R"irEdS,t) ^EdS rf"(fEds) ' 

where d{r) = p{r,t)+dp{t)~^/'^k{r), Jis a constant, p{t) is the EdS density, and p(r, t) = 3a{r)/8'KR^{r, t) is the initial 
averaged density. We see, therefore, that a general condition to have the constancy of R'" {r-Eds, t) / R" {^Eds, t) is that 
the second and third derivatives with respect to r of p{r,t) and kir) are proportional, so that the time dependence 
cancels. This is indeed the case in the AAG model. We remark that p{r,t) and k{r) can be seen as the defining 
functions for a LTB model, and this result holds thanks to the matching guaranteed by d'{r-Eds) = 0- 

2. Approximations for /3' 

We now discuss the approximation that af3' ~ R— ar. Starting from Eq. (|T7)) . with the approximation that W ~ 1, 
we have 



ri \ari a I 2 \ari a I X \ari 



1 /_R_ ir_ 

2 \ari a 



(A.7) 



Now from Fig. [21 the term in the square brackets oscillates about one, so it is not a bad approximation to set it to 
unity. The integral can then be done with the result that af3' ~ R— ar, used in Sect. lIVI 
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